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Abstract. - We studied phase separation in a particle interacting system under a large drive 
along X. We here identify the basic growth mechanisms, and demonstrate time self-similarity, 
finite-size scaling, as well as other interesting features of both the structure factor and the 
scaling function. We also show that, at late t in two dimensions, there is a unique t— dependent 
length increasing (.y (t) ~ t"-^'-" for macroscopic systems. Our results, which follow as a direct 
consequence of the underlying anisotropy, may characterize a class of nonequilibrium situations. 



Many binary mixtures, e.g. the alloy Al-Zn, which arc homogeneous at high temperature, 
develop coarsening macroscopic grains after a quench into the miscibility gap. The details of 
nucleation and spinodal decomposition as the system evolves towards coexistence of the two 
new phases determines various properties; e.g. hardness and resistivity of the alloy depend on 
how phase separation kinetics competes with the progress of solidification from the melt. 

The involved essential physics is rather well understood, partly due to computer simulation 
of lattice gases. [1, 2] An interesting task is now understanding more general situations, e.g. lack 
of isotropy, which bears great technological importance. Mixtures under a shear flow attracted 
considerable attention as a possible scenario. [3]- [5] We studied the kinetics of the driven lattice 
gas (DLG) [6] by extensive Monte Carlo (MC) simulation. This is appealing on several 
grounds. Firstly, the DLG is the most successful, microscopic metaphor for (nonequilibrium) 
anisotropic phenomena, and its time relaxation remains intriguing. [7]- [11] Furthermore, the 
common underlying anisotropy might induce essential macroscopic similarities between the 
DLG and the sheared lattice gas, [4] as recent results on critical behavior seem to suggest. [12, 11] 
In any case, the DLG is not affected by hydrodynamic effects, which makes physics simpler. [1] 

The DLG is a d— dimensional lattice gas (or, alternatively, binary alloy) at temperature T 
in which nearest-neighbor (NN) particle/hole exchanges are favored along one of the principal 
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lattice directions, say x.[7] There are variables = 1 (particle) or (hole) at each site 
z = 1, . . . , iV, a NN interaction according to H = —4 X^jvjv i^'i'^j ; toroidal boundary conditions, 
and a transition rate given (e.g.) by a —biased— Metropolis algorithm, uj{n n*) = 
mm{l,exp[— {AH + ES)/T]}, which conserves density p = N^^ J2i'>^i- Here, n* represents 
configuration n ={ni} after jumping of a particle to a NN hole. Ex may be interpreted as an 
(electric) field driving (charged) particles, AH = H{n*) — H{n), and 6 = (t1,0) for jumps 
along ±x or along any of the transverse directions, say y, respectively. 

The DLG was described as modelling fast ionic conduction, surface growth, traffic flow, 
etc. [7] A common feature of all these situations is anisotropy, and that steady states are out of 
equilibrium; both are essential features of the DLG induced by the rate uj. In general, uj violates 
detailed balance. This symmetry holds only for E — 0; the DLG then reduces to the familiar 
lattice gas with a unique (equilibrium) steady state independent of w(n — > n*). For any, even 
small E the steady state depends on w, and a different, qualitatively new behavior emerges. [12] 
As E is increased, one eventually reaches saturation (particles cannot jump backwards, — x), 
which is formally denoted as "E = oo" . 

For d = 2, p = I and E = oo (the only case to which we refer here for simplicity -also 
because this is a most interesting case [7]), a critical point occurs at T = ~ 1.4Tc {E = 0). 
Steady states below Tq' do not correspond to coexistence of two thermodynamic (equilibrium) 
phases as for E = 0. Instead, stable ordered states consist of one single stripe, corresponding to 
the liquid, rich-particle phase, and gas, poor-particle phase filling the remainder of the system. 
The interface is linear and flat, except for microscopic roughness (which slightly increases with 
decreasing E). One measures a net current of particles along x (its intensity increasing with 
T, and changing slope at T^).[7] 

Our simulations (fig-1) proceed by means of lo from a random state, until one or sometimes 
a few stripes are obtained. The code includes a list of X{t) particle-hole NN pairs from where 
the next move is drawn. Time is then increased by At = X{t)^^, so that its unit or MC step 
involves a visit to all sites on the average. Most evolutions are at T = 0.8Tc{E = 0) ~ 0.6T^ 
for which clustering is rather compact and, in practice, one can observe the full process of 
relaxation; we are assuming that, as observed for equilibrium states, [l]-[2] time evolution has 
the same properties in a wide region of the miscibility gap (see argument below) . The lattice 
is rectangular, x Ly, with sides ranging from 64 to 256 (and, exceptionally, 512). 

Starting from complete disorder, there is a very short regime in which small grains form (fig. 
la). Typical grains are anisotropic, stretched along x. One then observes a rapid coarsening 
to form macroscopic strings, fig. lb. We skip details concerning such nucleation and early 
phase separation. [13] We wish to notice, however, that sheared fluids seem to depict a similar 
behavior. That is, initial formation of small anisotropic clusters (inducing a larger growth rate 
along the flow than in the other directions), and stringlike domains extending macroscopically 
have been reported in fluids. [3, 5] 

After the initial (anisotropic) nucleation regime, strings coarsen further into well defined, 
relatively narrow stripes; fig.lc. The resulting multistripe states are not stable. They are 
partially segregated, and tend to relax towards the true stable state of one stripe. This may 
take a long macroscopic (say, infinite) time; in fact, the mean relaxation time increases with 
system size (see also [9]). As in the equilibrium case,[l] this is a consequence of the conservation 
of p, which makes the interfaces to depend on each other. It is true that certain individual 
runs sometimes block for a long time in a state with a few stripes; however, it does not seem to 
correspond to the average behavior. Typically, the number of stripes monotonically decreases 
with time, and the whole relaxation can easily be observed in computer simulations if one 
waits long enough (see caption for fig. 1). 

We shall assume that further coarsening occurs by monomer diffusion. Liquid stripes 
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(a) (b) (c) (d) 

Fig. 1. - A series of snapshots depicting growth at T ~ 0.6Tq' for Ly x Lx = 256 x 64 and t = 1 (a), 
50 (b), lO"* (c) and 2x10^ (d) MC steps. In tiiis particular run, two stripes were observed just before 
10^ MC steps while a single stripe was only reached after 10** MC steps. 

thus effectively diffuse and, eventually, collide and coalesce with one neighbor. This implies 
evaporation of a gas stripe. Therefore, given the particle/hole symmetry, our assumption is 
equivalent —though it allows for a more detailed description below— to assuming coarsening 
by stripe evaporation as in [9]. 

Let us evaluate sort of mobility Vi due to stripe diffusion via monomer events. Here Vi 
is the mean squared displacement of the stripe's center of mass per unit time. Consider a 
compact stripe of mean width C.{t) consisting of j = 1, . . . ,M particles at Ujit). Its center 
of mass is at VcmW = -^"^ Z]j yj(*)- Then Ve = Njne{{^ycraY) , with N^ne the frequency 
of events and ((AFcm)^) the mean squared displacement associated with one of them. [10] It 
ensues as the result of two competing processes: 

(A) A surface traps a monomer evaporated from the same interface. [14] Then Nme,A = 
exp(— 2/3Aj) where z/ is a frequency, the sum is over the surface particles, /3 is the 
inverse temperature and Aj is the number of broken bonds. For flat linear interfaces, Nme,A ~ 
'iuLx exp(— 2/?A) where A is the mean of Aj (we here multiplied by 2 to account for evaporation 
of surface holes reaching the surface again from the interior). One further has AYcm = M^^Sy, 
where Sy is the net displacement, and (for compact enough stripes) M x £{t), so that 

((An„)2) = {6y^) L-H-\ 
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Fig. 2. - The typical length I (t) versus t"^" for a lattice of size 64x64 (A), with = 3, and 256x64 
(O), with (j)y = ^. (Note that each set of data refers to different axis.) The inset shows Pearson Chi 
Square function, {4>y), for varying fits. 

(B) A hole jumps within the stripe, inducing AYf-m = 1/M or depending on the jump 
direction. One has Nme,B = '^'^Ph{T)L^iph{T), where phL^i is the number of holes in the 
stripe and ph is a jumping probabihty. At low T (small hole density ph), holes are isolated 
from each other, so that ph « 1. It ensues V^^^ 2iypfiL^^£^^ ■ 

As far as stripes are a distance £ apart from each other, they take a time = P' lT>i to 
meet, which increases width by I. Then dljdt ~ T>il~^ and, assuming that processes A and 
B are independent, 

Mjdt ~ L-'^{aAi~'^ + Q!B^"^) 

for low T and large E\ here, (xa = 4:i'{Sy'^)e~'^^^ and as = 2i/ph. 

This has some interesting implications. At late times, i{t) ~ at^^^ + aAj'^o-B with 
o? = SasL"^, i.e., hole diffusion (mechanism B) is dominant. Evaporation/condensation 
(mechanism A) results in behavior; this is predicted to matter earlier. The crossover 
between the two mechanisms is for t ~ Tcross X, I.e., a macroscopic, 

observable time. If we define the time at which a single stripe forms, £{tss) ~ ■^Ly, it 
follows that the t^^^ behavior is dominant for j{T, Ly, L^) = Tctoss/tss ^ 1- It also follows 
7^0 for finite T in the thermodynamic limit. Consequently, our theory predicts that the 
^1/3 growth is the general one to be observed, as it is also concluded in [9]; the t^/^ growth 
should only be observable in "small" —as defined below— systems. It is to be noted that, in 
equilibrium, the Lifshitz-Slyozov- Wagner mechanism (evaporation from small grains of high 
curvature and, after diffusion, condensation onto larger grains) implies i ~ t^^^, in accordance 
with experiments. 
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Fig. 3. - Scaling with size and time of the structure factor. The plot includes data from 64x64, 
128x128 and 256x256 lattices, and times t > 10"* MC steps. The relevant behaviors (see the main 
text) are indicated by dashed lines; k = kyi/L. 

Let us define a longitudinal length, ^ t'^=";[2, 9] one expects more rapid growth than 
transversely, (j)^ > 1/3. Concerning i— dependence, this length is only relevant initially, until 
well-defined stripes form. This interesting result concerns time dependence far from criticality. 
It is compatible with the possible existence of two correlation lengths near which describe 
thermal fluctuations. [12] Note also that the onset of the multistripe state, when only iy{t) is 
relevant, may be defined by Ixirms) = Lx'i i.e., Tyyig '"^ Lx 1 which is on the same macroscopic 
time scale as Tcross and Tss- 

The structure factor is, setting kx = (the dependence on hx is only relevant at early 
times): 



This soon develops a peak at ky = kmax{t) that shifts towards smaller wave number with 
increasing t. We define £s = 'i.Tr/kmax as a measure of the mean stripe width. Alternatively, 
one may use the first moment of S*, or the slope of the straight portion in a plot of h\[S{ky, t)] 
versus {ky — kmax)"^, which is known as Guinier radius . One may also monitor the number 
of stripes Ns and, averaging over all stripes in a given configuration, their maximum width, 
£max and £m = M/Lx- After performing an extra average over more than one thousand 
independent evolutions, all these quantities happen to exhibit essentially the same dependence 
on t; we denote t{t) this common behavior. 

To check our predictions, we plotted i{t) versus t'^'" for varying 4)y looking for the best 
linear fit. As illustrated in the inset of fig. 2, associated with this fit nicely confirms that 
0y 2± i for small systems {Ly < 128) while cjiy ^ for large ones [Ly > 256). This is further 
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confirmed by studying dln£ /dint for large £ (not shown). Our theory predicts that this 
crossover with size wiU occur for Tcross — '^ss' 

For flat linear interfaces, {Sy"^) ~ A < 6, 

Ph ~ exp(— 16/T) and ^ 1/q, the lattice coordination number. It then follows numerically 
that the crossover from the early (f>y = j to the general macroscopic | behavior will occur 
around Ly ~ 140 for = 64. Again, this is fully consistent with our observations. The 
reason behind is that surfaces (and, thus, mechanism A) dominate initially, more the smaller 
the system is; however, hole diffusion (B) tends to dominate as the system relaxes towards 
only two surfaces enclosing the whole of the liquid phase. 

The structure factor is an important tool for experimental analysis. Given that the 
DLG shows a unique f-dependent length after forming stripes, one should expect S{ky;t) oc 
£{t)F[ky£{t)]. This is confirmed in fig. 3. (A Ginzburg-Landau model for sheared mixtures 
has been shown to exhibit a similar property, though with two lengths both behaving dif- 
ferently than £{t) here. [5]) Fig. 3 also illustrates that data for different square lattices scale 
S£~^L~^ ^ F{riL~^), rj = ky£. A similar result does not hold for general, rectangular lattices, 
as one should probably expect due to more involved finite-size effects if Ly. 

The sphericallized function for a three-dimensional system relaxing to equilibrium was 
shown to satisfy S{k,t) = J{t) ■ F[kR{t)]. This turned out most useful given that R and 
J can simply be evaluated phenomenologically as the scaling lengths for the k and S axes, 
respectively. In particular, it then followed F{r]) = ^{r]) ■'^[r] ■ a{p,T)], with $ and ^ universal 
functions. $ describes the diffraction by a single grain, ^' is a grain interference function, and 
a characterizes the point in the phase diagram where the sample is quenched. In this way, 
it was shown [2] that » 1 except at small values of k, so that, for large r], F{r]) becomes 
independent of T, p, and even the substance investigated. Though we only obtained some 
weak, indirect evidence of this for the DLG (which was investigated systematically for just 
one phase point), it suggests that assuming a wide range of validity of our conclusions here is 
sensible. 

In equilibrium, $(ry) is predicted to decay according to Porod law,[l] ri~^ (for d = 2) at 
large rj. In the DLG, this concerns the region <§; ky^ <g; £{t), where ^y is the transverse 
correlation length, and the striped geometry implies a two-point correlation function C't{x, t) ^ 
|(1 — x£~^), X <^ £. Therefore, we predict that $ ~ ?7~^ is the characteristic anisotropic 
behavior corresponding to Porod law. This is confirmed in fig. 3. Interesting enough, this 
property of the scaling fmiction implies S ^ 1/ £ (t) ky for ky larger than for the Guinier 
gaussian region described above. Finally, we remark that the behavior $ ~ r]~^ is general 
except (for large systems) at large enough values of r] (fig. 3). In this case, the standard 
Porod law holds, which is induced by the very small, standard thermal clustering (which is 
here similar to the one in equilibrium). 

Summing up, we showed time self-similarity and simple finite-size scaling of the structure 
factor in a (nonequilibrium) particle model under a large drive following a quench to low T. In 
spite of some formal similarities, both the scaling function and the structure factor differ from 
the ones for the equilibrium (not driven) case. This is a consequence of the singular geometry 
(flat linear interfaces) induced by the drive. The scaling function $(ry) is secni to tend to an 
envelope ~ rj^+'^Z^ at small rj; it then follows a gaussian peak ~ exp[— const ■ {rj — rimaxy], 
and decays r]~^ and, finally, ry~^ at large r?. On the other hand, there is only one t— dependent 
relevant length, the mean interfaces distance. This generally grows as t^/^ due to hole diffusion 
in the bulk, though one may also observe t^^^ at very early times due to surface evaporation- 
condensation processes. The basic mechanisms are here much simpler than in equilibrium, 
again due to the underlying anisotropy. That is, the standard lattice gas typically develops 
grains of different sizes, showing a distribution from microscopic to macroscopic grains. This 
is the case for practically any t and, as a consequence, a detailed surface-tension-controlled 
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kinetic description is rather complex. [15] The DLG grain distribution, after a short transient 
time, exhibits a gap (at the, say, mesoscopic scale) between monomers and infinite grains, 
namely, stripes of varying width. Some of the DLG properties described here might hold in 
several other anisotropic situations. Wc hope related experiments will be performed that will 
motivate development of more complete theories. 
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